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Abstract 

Hybrid (exotic) mesons, which are important predictions of quantum chromodynamics 
(QCD), are states of quarks and anti-quarks bound by excited gluons. First principle lattice 
study of such states would help us understand the role of "dynamical" color in low energy 
QCD and provide valuable information for experimental search for these new particles. In this 
paper, we apply both improved gluon and quark actions to the hybrid mesons, which might be 
much more efficient than the previous works in reducing lattice spacing error and finite volume 
effect. Quenched simulations were done at f3 = 2.6 and on a £ = 3 anisotropic 12 3 x 36 lattice 
using our PC cluster. We obtain 2013 ± 26 ± 71 MeV for the mass of the 1 h hybrid meson 
qqg in the light quark sector, and 4369 ± 37 ± 99Mev in the charm quark sector; the mass 
splitting between the 1 h hybrid meson ccg in the charm quark sector and the spin averaged 
S-wave charmonium mass is estimated to be 1302 ± 37 ± 99 MeV. As a byproduct, we obtain 
1438 ± 32 ± 57 MeV for the mass of a P-wave 1++ uu or dd meson and 1499 ± 28 ± 65 MeV for 
the mass of a P-wave 1 ++ ss meson, which are comparable to their experimental value 1426 
MeV for the /i(1420) meson. The first error is statistical, and the second one is systematical. 
The mixing of the hybrid meson with a four quark state is also discussed. 
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1 Introduction 



The hybrid mesons qqg, which are bound states of quark q, anti-quark q and excited gluons g, 
have exotic quantum numbers inaccessible to qq mesons. There is increasing interest by BES and 
CLEO collaborations as well as experimentalists at Brookhaven and Jefferson labs in searching 
these new particles. Monte Carlo lattice technique is the ideal method not only for computing 
qq meson spectrum, but also for hybrids. Of course, the lattice QCD approach is not free of 
systematical errors. As is well known, the Wilson gauge and quark actions^ have 0(a 2 ) and 
O(a) errors respectively, where a is the lattice spacing. These errors are smaller only at very 
small bare coupling, and very large lattice volume is required to get rid of finite size effects. 

The idea of Symanzik improvement [2] is to reduce the lattice spacing errors, by adding new 
terms to the Wilson actions. Together with the tadpolejH] or non-perturbative improvement HJ , 
the pursuit of the Symanzik program has recently led to significant progress in lattice gauge 
theory, opening the possibility to approaching continuum physics on coarse and small lattices. 
Simulations on anisotropic lattices help getting very good signal in spectrum computations. 

There have been several quenched lattice calculations of hybrid meson masses, in either the 
light quark or heavy quark sector. In Refs.j3 EH d, the Wilson gluon action and quark action 
were used. In Refs.jSl the authors used Wilson gauge action and SW (clover) improved quark 
action |TT]. For the hybrid mesons containing heavy quarks QQg, Non-relativistic QCD (NRQCD) 
quark action [TUl I12j and LBO (leading Born-Oppenheimer) quark action t f3] were also applied to 
the study of charmonium. As mentioned in B,efs.[T H ITn l ITH]. these approaches have advantages 
and disadvantages. 

In this work, both improved gluon and quark actions on the anisotropic lattice are employed, 
which should have smaller systematical errors, and should be more efficient in reducing the lattice 
spacing and finite volume effects. We will present results for the 1 ^ hybrid mass and the 
splitting between the 1 ^ hybrid mass and the spin averaged S-wave mass for charmonium. As a 
byproduct, we will also estimate the /i(1420) mass. 

The remaining part of the paper is organized as follows. We describe the improved gluon and 
quark actions in SecEl and operators for the qq mesons and qqg mesons in Sec|21 Simulation 
details are given in Sec0]and results and conclusions are summarized in SecJSl 

2 Improved actions on the anisotropic lattice 

The total lattice action is S = S g + S q . S g is the classical-improved gluon action [T8| IT9l 120] : 

x,j<k 3 x,j 

which differs from the continuum Yang-Mills action in imaginary time J d A x F™F"/4 by 0(a 2 ) 
and O(aj) at classical level, with at and a s being the temporal and spatial lattice spacing respec- 
tively, g is the gauge coupling; x is the lattice site and the subindeces j, k and 4 in Eq. Q are 
respectively the spatial and temporal directions; £o is the bare anisotropic parameter. P stands 
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for a 1 x 1 plaquette and R for a 2 x 1 rectangle. More explicitly, they are 
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where Uj(x) is the gauge link variable at site x and j-th direction. The clover action for quarks 
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which differs at classical level from the continuum Dirac action J d 4 x q^^D^ + m q )q by 0(a 2 ). 
(There is another way to reduce the lattice spacing error using the Hamber-Wu quark action [T7|.^ 
The continuum quark field q is related to the lattice counter part ip by tp = y / a^/(2K t )q, and the 
bare quark mass is related to Kt and k s by 

a t m q = l/(2«t) - 3r s K s /K t - 1. (4) 

For the gauge field tensor JP we use the clover deaf construction |llU21| . and the clover coefficients 
at classical level are 

1 + Cor s 



C s = r s , C t 



(5) 



where r s is the Wilson parameter and it is usually taken to be 1. 

The tadpole-improved actions are more continuumdike, and can be obtained by following 
replacement 
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£ is just the aspect ratio of spatial and temporal lattice spacings a s /at\ u s and ut are tadpole 
parameters are taken to be 
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Table 1: Source and sink operators used in Eq.JSJ), the two point Green function of mesons. 

3 Meson operators and green functions 

The hybrid operators are constructed by combining a quark, an anti-quark and the chromo- 
electric or chromo-magnetic field to form a color singlet with the given spin J, parity P and 
charge conjugation C. The generic structure is i[) ( j 1 Ttp'j 2 2 F ClC ' 2 , where c\ and C2 stand for color 
indices, and f\ and f% for flavor indices; V is some combination of Dirac matrices and spatial 
derivatives. The details are given in Ref.:7j. For convenience, we present in Tab. ^ various source 
and sink operators for qq mesons, qqg mesons and a four quark state we have used. We consider 
such a q 4 operator (four quark state) because it can mix with the hybrid mesons through the 
hairpin diagram. 

In TabQ the chromo-electric and chromo-magnetic fields bring in J PC = 1 and 1^ re- 
spectively, eijk'ip^j dk ift is a P-wave operator for f\ with J PC = 1 ++ . The operator di=di — di 
have J PC = 1 , where C = — 1 is because that C interchanges the quark and the anti-quark. 

To use information from the different source and sink operators, we fit simultaneously the 
correlation function of two or more operators with different amplitudes A{ and Aj for each source 
Oi and sink Oj listed in Tab^ but with the same masses m e g for all of them. For large t, it 
behaves as 

J2 < 0,(6, 0)O j(x,t) >= AiAje-™**** + ... (8) 

X 

4 Simulation results 

Part of our work is based on the MILC public code 22 , which was written for isotropic lattices. We 
modified the code to anisotropic lattices. To check the new code for clover quarks, we computed 
some charmonium masses using the gauge configurations generated by the standard one-plaquette 
gauge action, and compared with those of the CP-PACS collaboration 16 . To check the new code 
for gluons with the tadpole improved gauge action, we calculated the tadpole parameter u s and 
compared results with those of Alford et al. 2Qj. The simulations were performed on our PC cluster 
for parallel computing [221 122 ESI EH1 Elj The results are well consistent with the literature, and 
details of comparison will be reported elsewhere. 
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Table 2: Simulation parameters for hybrid meson spectrum calculations. 

The results presented in this paper use the SU(3) pure gauge configurations generated with 
the tadpole-improved gluon action 



R l V" f 5P J,k 1 Rj,h 1 Rk,j\ at\^( 4P 3A 1 R iA (q\ 
V t I Q „,4 1 9 „,6 10 J ^Z^ U „2 n „4 ' W 



3 ^£ 4^, V3 it* 12 u% 12 «6 y ^ S ^V3 12 «4 

using Cabibbo-Marinari pseudo-heatbath algorithm. The configurations are decorrelated by SU(2) 
sub-group over-relaxations. We calculated the tadpole parameter u s in a self-consistent manner. 
At j3 = 2.6 and £ = 3, we generated 90 pure gauge configurations on the 12 3 x 36 lattice. Although 
such an ensemble is not very big, it is bigger than earlier simulations by UKQCD and MILC 
collaborations [HI 13 IH1 El on isotropic lattices. To practically implement the tadpole-improvement 
for the quarks, it is more convenient to tune the parameters Kt, k s (or ( = Kt/ K s ), C^ 1 and Cj l 
of the following quark action 



S q = $(?)if>(x) - n t - 74)*7 4 (^(x + 4) + ${x)(l + lA )Ut (x)if>(x - 4)] 

X X 
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+ iKgCj 1 Y ip(x)a jk F jk (x)ip(x) +iK s C^ l Y%l)(x)a jA F ji {x)tp{x), (10) 

x,j<k x,j 

where the clover constants are 

without rescaling the link variables. The simulation parameters are listed in Tab|2j 

We obtained the quark propagator by inverting the quark matrix A in the tadpole-improved 
Sq = ^2x,y^( x )^x,y^{y) by means of BICGStab algorithm. The residue is of O(10 -7 ). We 
evaluated the correlation functions with various sources and sinks listed in Tab. ^ at four values 
of the Wilson hopping parameter (kj = 0.4119, 0.4199, 0.4279, 0.4359). 

The effective mass is obtained by fitting the correlation function to Eq.JSJ). In Fig^ we show 
the effective masses for the tt, p, fi ordinary mesons and 1 ^ exotic meson at Kt = 0.4359. For 
the ordinary mesons, we used the their corresponding operators in Tab. ^ as both source and 
sink. For the exotic meson, we tried two different cases: (1) the 1 h operator as both source and 
sink; (2) the q 4 source and 1 h sink, which give consistent results within error bars. 

In Tab |3J we present the data for qq meson masses at four different Kt, and in Tab|I] the 
data for the 1 ^ exotic meson. Fig|2] plots the pion mass squared as a function of l/re*. By 
linearly extrapolating the data to the massless limit, we obtain ^ hiral = 0.5945(3)(11), which 
is also shown in the figure. The regression functionality of the software Xmgr is used in the 
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0.727(7) 
0.711(9) 
0.694(8) 
0.679(5) 
0.458(3)(15) 


3.5/7 
7.8/7 
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fi P-wave 
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0.4279 
0.4359 
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3-11 
3-11 
3-11 
3-11 
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1.069(17) 
1.052(13) 
1.043(14) 
1.030(9) 
0.857(14)(27) 


4.9/5 
4.6/5 
8.9/5 
1.2/5 



Table 3: Effective masses of ordinary qq mesons at different Kt- 

extrapolation. In FigEl we show the effective mass of other mesons as a function of 1/Kt an d 
extrapolation to the chiral limit. In the chiral limit, we have a t m p (n t -» «f iral ) = 0.458(3)(15). 
Using the experimental value m p = 768.5MeV, we get of 1 = 1.679(2) (7)GeV, and 

a t O.llFermi, a s « 0.33Fermi. (12) 

The UKQCD, MILC and CP-PACS collaborations used the unimproved Wilson gauge action 
to generate configurations. They had to work on very large /?(> 6), corresponding to very small 
a s {< 0.1 Fermi), to get rid of the finite spacing errors. They had also to use very large lattices 
L 3 > 20 3 , to avoid strong lattice size effects at such small a s . In comparison, our lattices are much 
coarser, and the number of lattice sites is much smaller. (Of course, finite size effects could be 
ignored because the physical size of our spatial lattice 12 3 a 3 = 3.96 3 Fermi 3 should be big enough.) 
Our results for the effective mass indicate the existence of a very wide plateau, which are better 
at least than the previous work on isotropic lattices. 

5 Predictions of masses and discussions 

In this concluding section, we predict the 1 *~ hybrid meson masses and discuss briefly the sys- 
tematical errors. 

By extrapolating the effective mass of the 1 ^ hybrid meson (TabHJ) to the chiral limit, and 
using at determined in Sec. |1J we get 2013 ± 26 ± 71 Mev. In Tab. [SJ we compare the results from 
various lattice methods. Our result is consistent with the data of MILC collaboration who 
obtained using the Wilson gluon action and clover quark action on much larger isotropic lattices 
and much smaller a s . 
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Source(s)— >Sink 


Fit Range 


X 2 /dof 


atm 


0.4119 


1-+ -> 1-+ 


3-10 
3-10 


6.3/5 
3.6/5 


1.466(26) 
1.459(25) 


0.4199 




3-10 
3-10 


6.8/5 
5.5/5 


1.456(28) 
1.449(36) 


0.4279 


1-+ -> 1"+ 


3-10 
3-10 


1.0/5 
3.8/5 


1.437(28) 
1.434(31) 


0.4359 


1-+ -> 1-+ 

<fM-+-i-+ 


3-10 
3-10 


3.3/5 
7.0/5 


1.418(35) 
1.400(37) 


K cbrial = 0.5945 

K charm _ Qjggg 


1-+ 1-+ 
1-+ -► 1-+ 


Extrapolation to the chiral limit 
Tuning 




1.197(8)(26) 
2.603(12)(23) 



Table 4: Effective mass of the exotic 1 + meson at different Kf. 



Also in Tab. we present our results for the 1 + hybrid meson mass in the charm quark 
sector, using the method discussed in Refs. [7j|3]. The corresponding K^ harm = 0.1806(5)(18) is 
obtained by tuning {m n (K t -> 4 harm ) + 3m p (K t -> 4 harm ))/4 = (rn Vc + 3m J/ , / ,)/4=3067.6MeV, 
where on the right hand side, the experimental data m Vc = 2979.8 Mev and mj/^ = 3096.9 MeV 
are used. The 1 ^ hybrid meson mass at our 1/K£ harm is m\-+ = 4369 ± 37 ± 99 MeV, which 
is consistent with the MILC datajjj. The splitting between the hybrid meson mass and the spin 
averaged S-wave mass [mx-+ — (m Vc + 3mj/^)/4], at our K £ harm i s 1302±37±99 MeV, consistent 
with the CP-PACS data, obtained using the Wilson gluon action and NRQCD quark action on 
much larger anisotropic lattices and much smaller a s . 

As a byproduct, in Tab. El we give the f\ P-wave 1 ++ meson in the chiral limit, as well as 
their experimental values [25]. If we assume that the pion is massive and /i(1420) is made of ss, 
the fi P-wave 1 ++ meson mass would be 1499 ± 28 ± 65 MeV. 

In conclusion, we have used the tadpole-improved gluon action and clover action to compute 
the hybrid meson masses on much coarse anisotropic lattices. The main results are summarized 
in Tab|S]and compared with other lattice approaches. In our opinion, our approach is much more 
efficient in reducing systematic errors due to finite lattice spacing and finite volume. 

In our study, the quoted errors are statistical and systematical. We have not considered the 
quenching errors. (Only in the case of unimproved Wilson gluon action plus the clover quark 
action, the effects of dynamical quarks have been studied 29 a by applying the algorithm of Refs. 
[^ISOl, and are found to have little influence on the spectrum of bottomonium spectrum.) We 
have not taken into account the contamination of excited states when fitting the correlation 
function; this is a common problem in lattice Lagrangian formulation; Monte Carlo Hamiltonian 
method j3T | 132* 1 133] proposed recently might be helpful in this aspect. Also, we have assumed that 
with the improved gluon and quark actions, the dependence of the results on the lattice spacings 
is negligible, and we haven't extrapolated the results to the continuum limit. We hope to discuss 
these issues in the future. 
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Light l~ + qqg (GeV) 


Method 


Ref. 


1.97(9)(30) 
1.87(20) 
2.11(10) 
2.013(26)(71) 


Isotropic S g (W) + S q {W) 
Isotropic Sj^(W) + S^(SW) 
Isotropic Sj^W) + SJ T (SW) 
Anisotropic Sj l (l x 1 + 2 x 1) + Sj ! (SW) 


MILC970 
UKQCD97[H] 
MILC990 
ZSU (this work) 


\- + ccg (GeV) 


Method 


Ref. 


4.390 (80) (200) 
4.369 (37) (99) 


Isotropic S g (W) + S q (W) 
Anisotropic Sj l (l x 1 + 2 x 1) + Sj 1 (SW) 


MILC97[Z1 
ZSU (this work) 


1 v ccg -IS cc splitting (GeV) 


Method 


Ref. 


1.34(8)(20) 
1.22(15) 
1.323(13) 
1.19 
1.302(37)(99) 


Isotropic S g (W) + S q {W) 
Isotropic Sj : (W) + 5j!(SW) 
Anisotropic sf(W) + ^(NRQCD) 
Isotropic Sj\l x 1 + 2 x 1) + 5^(LBO) 
Anisotropic Sj l (l x 1 + 2 x 1) + Sf^ (SW) 


MILC97[Zj 
MILC990 
CP-PACS99QU1 
JKM99 [E] 
ZSU (this work) 



Table 5: Predictions for the masses of 1 ^ hybrid mesons qqg in the light quark sector and ccg in 
the charm quark sector. Different lattice approaches are compared. Abbreviations: W for Wilson, 
1x1 + 2x1 for the plaquette terms plus the rectangle terms, SW for Sheikholeslami-Wohlert 
(Clover), TI for tadpole-improved, NRQCD for non-relativistic QCD, and LBO for leading Born- 
Oppenheimer. 



/i P-wave (GeV) 


Method 


Ref. 


Experiment 


1.438 (32) (57) 


Anisotropic S£' x (l x 1 + 2 x 1) + S'^(SW) 


ZSU (this work) 


1.426 /i(1420)? 



Table 6: Predictions for the masses of the 1 ++ P-wave mesons qq. 
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Figure 1: The effective masses for the it (triangle down), p (circles), f\ P-wave (square) mesons 
and 1 h exotic meson (the diamond for 1~ + source and the triangle up for the q A source) at 
K t = 0.4359. The error bars reflect the sum of statistical and systematic errors. 




Figure 2: Pion mass squared as a function of 1/nf The error stands for the sum of statistical 
and systematic errors. 
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Figure 3: Effective masses as a function of and extrapolations to the chiral limit. The circles 
are p masses for both p source and sink, the squares are masses for the f\ P-wave operator as 
both source and sink, and the triangles are masses for the 1 _+ — > 1 h hybrid operator. 
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